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ABSTRACT 

With the aim of determining if Milky Way (MW) progenitors could be identified as 
high redshift Lyman Alpha Emitters (LAEs) we have derived the intrinsic properties 
of z sa 5.7 MW progenitors, which are then used to compute their observed Lya lumi- 
nosity, La, and equivalent width, EW. MW progenitors visible as LAEs are selected 
according to the canonical observational criterion. La > iC^ erg s~^ and EW > 20 
A. Progenitors of MW-like galaxies have La = lO'^^"^'^-^^ erg s~^, making some of 
them visible as LAEs. In any single MW merger tree realization, typically only 1 
(out of « 50) progenitor meets the LAE selection criterion, but the probability to 
have at least one LAE is very high, P = 68%. The identified LAE stars have ages, 
t^, ~ 150 — 400 Myr at z « 5.7 with the exception of five small progenitors with 
i* < 5 Myr and large EW = 60 - 130 A. LAE MW progenitors provide > 10% of the 
halo very metal-poor stars [Fe/H]< —2, thus establishing a potentially fruitful link 
between high-z galaxies and the Local Universe. 

Key words: cosmology: theory - galaxies: high redshift, evolution, stellar content, 
luminosity function - stars: formation, population II - 



1 INTRODUCTION 

Lyman Alpha Emitters (LAEs) are galaxies identified by 
means of a very strong Lya line (1216 A). Advances 
in instrument sensitivity and specific spectral signatures 
(strength, width and continuum break bluewards of the line) 
have enabled the confirmation of hundreds of LAEs in a 
wide redshift range, at z ~ 2.25 (Nilsson et al. 2008), z « 3 
(Cowie & Hu 1998; Steidel et al. 2000; Matsuda et al. 2005; 
Venemans et al. 2007; Ouchi et al. 2008), z « 4.5 (Finkel- 
stein et al. 2007), z « 5.7 (Malhotra et al. 2005; Shimasaku 
et al. 2006) and z ~ 6.6 (Taniguchi et al. 2005; Kashikawa 
set al. 2006). 

LAEs have by now been used extensively as probes of 
both the ionization state of the intergalactic medium (IGM) 
and probes of high redshift galaxy evolution (Santos 2004; 
Dayal, Ferrara & Gallerani 2008; Nagamine et al. 2008; 
Dayal et al. 2009; Dayal, Ferrara & Saro 2010; Dayal, Maselli 
& Ferrara 2010). However, in spite of the growing data sets, 
there has been no effort to establish a link between the prop- 
erties of these early galaxies to observations of the local Uni- 
verse, in primis the Milky Way (MW). Our aim in this work 
is to investigate the possible connection between the Galac- 
tic building blocks and LAEs at a time when the Universe 
was « 1 Gyr old. This will allow us to answer to questions 
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as: are the progenitors of MW-like galaxies visible as LAEs 
at high redshifts? How can we discriminate amongst LAEs 
which are possible MW progenitors and those that are not? 
What are the physical properties of these Galactic building 
blocks? 

To this end, we adopt a novel approach of coupling the 
semi-analytical code GAMETE (Salvadori, Schneider & Fer- 
rara 2007; Salvadori, Ferrara & Schneider 2008, Salvadori & 
Ferrara 2009), which traces the hierarchical build-up of the 
Galaxy, successfully reproducing most of the observed MW 
and dwarf satellite properties at z = 0, to a previously de- 
veloped LAE model (Dayal, Ferrara & Gallerani 2008; Dayal 
et al. 2009; Dayal, Ferrara & Saro 2010), that reproduces a 
number of important observational data sets accumulated 
for high-z LAEs. 



2 OBTAINING THE MW PROGENITORS 

We start by summarizing the main features of GAMETE, which 
is used to build-up 80 possible hierarchical merger histo- 
ries and to derive the properties of the MW progenitors at 
z = 5.7' . First, the possible hierarchical merger histories of 



T We specifically choose z = 5.7 for all our calculations as it 
represents the highest redshift for which a statistically significant 
sample of confirmed LAEs is available. 
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Figure 1. Cumulative Lyo luminosity function at redshift z fs 
5.7. The solid line shows the mean Lya LF of the MW progeni- 
tors averaged over 80 realizations of the merger tree; the dashed 
area represents the ztlfr spread among different realizations. For 
comparison only we show (points) the field LAE LF (Shimasaku 
et al. 2006). 



a MW-size dark matter (DM) halo are reconstructed up to 
z — 20 via a Monte Carlo algorithm based on the extended 
Press-Schechter theory (see Salvadori, Schneider & Ferrara 
2007 for more details) . The evolution of gas and stars is then 
followed along each hierarchical tree by assuming that: (a) 
the initial gas content of DM haloes is equal to the universal 
cosmological value [Qi,/^m.)Mh, vyhere M^ is the DM halo 
mass, Qb {^Im.) is the baryonic (DM) density parameter; (b) 
at any redshift, there exists a minimum halo mass to form 
stars, Msf{z), whose evolution accounts for the suppression 
of star formation (SF) in progressively more massive objects 
due to radiative feedback effects (see Fig. 1 of Salvadori & 
Ferrara 2009); (c) the gradual accretion of cold gas, Mc, into 
newly virializing haloes is regulated by a numerically cali- 
brated infall rate (Keres et al. 2005); (d) the SF rate, M,, 
is proportional to the mass of cold gas inside each galaxy, 
M* — etMc/tff, where e, is the SF efficiency and f// the 
halo free fall time; (e) in halos with a virial temperature, 
T^ir < 10* K (minihalos), the SF efficiency is reduced as 
e = e4l + {T^ir/2 X 10*K)"^]"^ due to ineffective cooling by 
H2 molecules. The chemical enrichment of gas, both in the 
proto-Galactic halos and in the MW environment is followed 
simultaneously by taking into account the mass-dependent 
stellar evolutionary timescales and the effects of mechani- 
cal feedback due to supernova (SN) energy deposition (see 
Salvadori, Ferrara & Schneider 2008 for more details). 

The two free parameters of the model (star formation 
and wind efficiencies) are calibrated by reproducing the 
global properties of the MW (stellar/gas mass and metal- 
licity) and the Metallicity Distribution Function (MDF) of 
Galactic halo stars (Salvadori, Schneider & Ferrara 2007, 
Salvadori, Ferrara & Schneider 2008); Msf{z) is fixed by 
matching the observed iron-luminosity relation for dwarf 



spheroidal galaxies (Salvadori & Ferrara 2009). They are as- 
sumed to be the same for all progenitors in the hierarchical 
tree. 



3 IDENTIFYING LAES 

By using GAMETE we obtain the total halo/stellar/gas masses 
{Mh, Mf, Mg), the instantaneous SF rate (M,), the mass 
weighted stellar metallicity (Z,) and the mass- weighted stel- 
lar age (i,) of each MW progenitor, in each of the 80 real- 
izations considered. These outputs are used to calculate the 
total intrinsic Lya (I/™*) and continuum luminosity (L™*) 
which include both the contribution from stellar sources and 
from the cooling of collisionally excited neutral hydrogen 
(H I ) in the interstellar medium (ISM) (Dayal, Ferrara & 
Saro 2010). 

The intrinsic Lya luminosity can be translated into the 
observed luminosity such that La ~ U^^ faTa, while the ob- 
served continuum luminosity, Lc is expressed as La = L^T'^fc- 
Here, fa (fc) are the fractions of Lya (continuum) photons 
escaping the galaxy, undamped by the ISM dust and Ta 
is the fraction of the Lya luminosity that is transmitted 
through the IGM, undamped by H i 1^. 

The main features of the model used to calculate /c, fa 
and Ta are; (a) for each MW progenitor the dust enrich- 
ment is derived by using its intrinsic properties (M*, t., Mg) 
and assuming Type II supernovae (SNII) to be the primary 
dust factories. The dust mass, Md, is calculated including 
dust production due to SNII (each SN produces O.5A/0 of 
dust), dust destruction with an efficiency of about 40% in 
the region shocked to speeds > 100 kms~^ by SNII shocks, 
assimilation of a homogeneous mixture of dust and gas into 
subsequent SF (astration), and ejection of a homogeneous 
mixture of gas and dust from the galaxy due to SNII. (b) 
Following the calculations presented in Dayal, Maselli & Fer- 
rara (2010) for high-z LAEs, the dust distribution radius, ra, 
is taken to scale with the gas distribution scale, rg, such that 
ra ~ 0.5rg ■^ (c) fc is calculated assuming a slab-like dust 
distribution and we use fa = 1.3/c, as inferred for LAEs 
at z f» 6. (d) Ta is calculated using the mean photoioniza- 
tion rate predicted by the Early Reionization Model (ERM, 
reionization ends at z; ~ 7) of Gallerani et al. (2008), accord- 
ing to which the neutral hydrogen fraction xhi = 7 x 10~^ 
at 2 ~ 5.7. Complete details of these calculations can be 
found in Dayal, Ferrara & Saro (2010) and Dayal, Maselli & 
Ferrara (2010). 

Progenitors are then identified as LAEs based on the 
currently used observational criterion; La > W^^ ergs~^ and 
the observed equivalent width, La/Lc > 20 A. We use a 
comoving volume of the Milky Way environment (SOMpc'^) 
to calculate the number density of the progenitors visible 
as LAEs and the observed LAE luminosity function (LF), 
shown in Fig. 1. 



■f The continuum band (1250-1500 A) is chosen so as to be unaf- 
fected by H I absorption. 

3 the gas distribution radius is calculated as rg = 4.5Ar200i where 
the spin parameter, A = 0.05 (Ferrara, Pettini & Shchekinov 
2000) and r200 is the virial radius. 
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Figure 2. Probability of finding a number, Nla.Es, of LAEs in 
a single MW realization, averaged over 80 hierarchical merger 
histories at z Ri 5.7. 



4 RESULTS 

We start by comparing the number density of field LAEs 
observed at z ~ 5.7 with the average LF of the Milky Way 
progenitors at the same epoch (Fig. 1). The number den- 
sity of MW progenitors decreases with increasing Lya lu- 
minosity, reflecting the higher abundance of the least mas- 
sive/luminous objects in ACDM models. The MW progen- 
itors cover the entire range of observed Lya luminosities. 
La = 10*^~*"^'^^erg s~^. We then conclude that among the 
LAEs observed at z ^ 5.7 there are progenitors of MW-like 
galaxies. For any given La the number density of MW pro- 
genitors is higher than the observed value because the MW 
environment is a high-density, biased region. As discusses in 
the final Section, uncertainties on the treatment of dust may 
also play a role. 

In Fig. 2 we show the probability distribution func- 
tion (PDF), P, of finding a number Nlaes of LAEs in 
any given MW realization. The PDF has a maximum (P — 
0.42) at Nlaes ~ 1; while in any single MW realization, 
there are ~ 50 star-forming progenitors with stellar masses 
M* ^ lO' Mq, on average only one of them would be visible 
as a LAE. We note that P = 0.32 for Nlaes = 0, while 
it rapidly declines (P < 0.2) for Nlaes > 1. We conclude 
that the MW progenitors that would be observable as LAEs 
Sit z K, 5.7 are rare (« 1/50), but the probability to have 
at least one LAE in any MW hierarchical merger history is 
very high, P = 68%. 

Let us now consider the physical properties of the build- 
ing blocks of the MW. M* represents the dominant physical 
factor to determine whether a progenitor would be visible as 
a LAE, since it governs the intrinsic Lya/continuum lumi- 
nosity, the dust enrichment (and hence absorption) and Ta, 
as both the size of the ionized region around each source and 
the H I ionization fraction inside it scale with M« (Dayal, 
Ferrara & Gallerani 2008). This implies the existence of a 



SF rate threshold for MW progenitors to be visible as LAEs 
which is Mr'" « O.gMeyr-i (panel (a) of Fig. 3). Since 
M« oc Mg (see Sec. 2), such a lower limit can be translated 
into a threshold gas mass: M™'" « 8x lO'^M© (panel (b)). In 
turn, the gas content of a (proto-) galaxy is predominantly 
determined by the assembling history of its halo and the ef- 
fects of SN (mechanical) feedback. While the most massive 
MW progenitors [Mh > 10^° Mq) display a tight Mg - Mh 
correlation, the least massive ones are highly scattered, re- 
flecting the strong dispersion in the formation epoch/history 
of recently assembled halos (panel (b)). As a consequence, 
LAEs typically correspond to the most massive progenitors 
of the hierarchical tree, i.e. the major branches (black points 
in the panels). In particular, we find that all haloes with 
Mh > 1O^°M0 (40 haloes) are LAEs. At decreasing Mh, in- 
stead, the progenitors can be visible as LAEs only by virtue 
of a high gas mass content or extremely young ages; there 
are 5 such objects, with Mh ~ 10^ Mq, as seen from panel 
(b). 

By comparing the panels (b) and (c) of Fig. 3, we can 
see that Md of the Galactic building blocks closely tracks 
Mg. Since the gas mass content of the z « 5.7 MW pro- 
genitors is reduced by ~ 1 order of magnitude with respect 
to the initial cosmic value, due to gas (and dust) loss in 
galactic winds, their resulting dust mass is relatively small: 
Md ~ 10''"^'^ M0. As a consequence all the progenitor galax- 
ies have a color excess E(B-V)< 0.025. In panel (d) we see 
that as expected, E(B-V) increases with M«, i.e. more mas- 
sive galaxies are redder. However the trend is inverted for 
M* < 10* M©; even though the dust masses in these low 
mass objects {M < 10^ Mq) is Md ^ W^'^Mq, due to their 
small virial radius, both the gas distribution scale, and hence 
the dust distribution scale are very small (see Sec. 3). The 
concentration of the dust in a small area leads to a large dust 
attenuation and hence, a large value of the color excess. 

MW progenitors visible as LAEs are generally interme- 
diate age objects, i, « 150 — 400 Myr, as seen from panel 
(e); the largest progenitors tend to be the oldest ones, an ex- 
pected feature of standard hierarchical structure formation 
scenarios. However, the ages show a large scatter, especially 
at decreasing M« , reflecting the great variety of assembling 
(and SF) histories of recently formed halos. Interestingly the 
five Mh ~ 10^ M© newly formed {z < 6) progenitors visible 
as LAEs have a very young stellar population, t« < 5 Myr. 
The high AL induced by the large mass reservoir and the 
copious Lya production from these young stars makes them 
detectable. 

In panel (f) we can see that these newly virializing 
galaxies are metal-poor objects, with an average stellar 
metallicity Z ~ 0.016 — 0.044Zq. As these galaxies host 
a single and extremely young stellar population, such low 
Z value reflects the metallicity of the MW environment at 
their formation epoch, z « 5.7 (see the middle panel in Fig. 1 
of Salvadori, Ferrara & Schneider 2008). The more massive 
MW progenitors visible as LAEs, instead, are more metal 
rich, ^ ~ 0.3 — IZq; their intermediate stellar populations 
form during a long period (see Fig. 4) and from a gas that 
is progressively enriched by different stellar generations. 

The fact that the physical properties of the LAEs pro- 
genitors obtained with our model are consistent with those 
inferred from the observations of field LAEs (Fig. 3) is a 
notable success of our model. 
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Figure 3. Physical properties of z fs 5.7 MW progenitors in 80 different hierarchical merger histories. We show: (i) all the progenitors 
(yellow circles), (ii) the major branches of each hierarchical tree (black filled points) and (iii) the progenitors identified as LAEs (colored 
open symbols). LAEs pertaining to the same (different) realizations are shown with the identical (different) colored symbols (see the 
text and Fig. 2). Triangles are used for those realizations in which there is only one LAE. As a function of the total stellar mass Af* 
the various panels show: (a) the instantaneous star formation rate, Af«; (c) the dust mass, M^; (d) the color excess E{B — V); (e) the 
average stellar age, i,; (f) the average stellar metallicity Z. Panel (b) shows the relation between the halo and gas mass, with the cosmic 
value {Qi,/Qrn)Mh pointed out by the dashed line. Points with errorbars are the observational LAE data collected by Ono et al. 2010 
(magenta squares, 1 LAE, four different models), Pirzkal et al. 2007 (blue triangles, 3 LAEs) and Lai et al. 2007 (green circles, 3 LAEs). 
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Figure 4. Star formation history of a typical MW progenitor 
identified as a LAE at 2 ss 5.7. 
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Figure 5. Fraction of LAE MW halo relic stars as a function of 
their iron abundance. The histogram shows the average fraction 
among 80 realizations of the hierarchical tree; the shaded area 
shows the ±lcr dispersion among different realizations. 



As mentioned above, the scatter in the LAE ages orig- 
inates from different assembling and SF histories (SFH) of 
the MW progenitors. This can be better understood by con- 
sidering the SFH of a typical '' MW progenitor identified as a 
LAE, as shown in Fig. 4. We find that since most LAEs typi- 
cally correspond to the major branch, their progenitor seeds 



il We define as "typical" a LAE whose properties match the av- 
erage values: M^, ^ 10^°Mq, Mg ?d 10* Mq, M* ^ 2.3 Mq/jt, 
t, ^ 230 Myr. 



are associated with high-o- peaks of the density field virial- 
izing and starting to form stars at high redshifts (z ~ 16). 
The SF rapidly changes in time, exhibiting several bursts 
of different intensities and durations, which follow merging 
events refueling gas for SF. During the "quiet" phase of ac- 
cretion, instead, SN feedback regulates SF into a more gentle 
regime. The duration and intensity of the peaks depends on 
the effectiveness of these two competitive physical processes. 
At high redshifts {z > 10) the peaks are more pronounced 
because (i) the frequency of major merging is higher and 
(ii) mechanical feedback is stronger given the shallower po- 
tential well of the hosting halos (M ~ 10*'^ Mq, Salvador!, 
Ferrara & Schneider 2008). 

We finally turn to the last question concerning the con- 
tribution of old and massive progenitors seen as LAEs, to 
the very metal-poor stars ([Fe/H]< —2) observed in the MW 
halo. In Fig. 5 we show the fractional contribution of long- 
living stars from LAEs to the MDF at z; = 0. LAEs provide 
> 10% of the very metal-poor stars; the more massive the 
LAE, the higher is the number of [Fe/H] < — 2 stars they con- 
tribute. This is because such metal-poor stellar fossils form 
at 2 > 6 in newly virializing halos accreting pre-enriched 
gas out of the MW environment (see Salvador!, Schneider & 
Ferrara 2007 and Salvador! et al. 2010). By z « 5.7, many of 
these premature building blocks have merged into the major 
branch, i.e. the LAE. Because of the gradual enrichment of 
the MW environment, which reaches [Fe/H]~ —2 at z « 5.7 
(see Fig. 1 of Salvador!, Ferrara & Schneider 2008), most of 
[Fe/H]> —2 stars form at lower redshifts, z < 5.7, thus pro- 
ducing the drop at [Fe/H]> —2. Note also the rapid grow of 
^LAEa/^MW ^^ [Fe/H] > —1, which is a consequence of the 
self-enrichment of building blocks resulting from internal SN 
explosions. 



5 CONCLUSIONS 

We have linked the properties of high-jz LAEs to the Local 
Universe by coupling the semi-analytical code GAMETE to a 
previously developed LAE model. According to our results 
the progenitors of MW-like galaxies cover a wide range of 
observed Lya luminosity, La = IQ^^"*'^-^^ erg s~^, with La 
increasing with M, (or, equivalently, Mh); hence some of 
them can be observed as LAEs. In each hierarchical merger 
history we find that, on average, only one star-forming pro- 
genitor (among ~ 50) is a LAE, usually corresponding to 
the major branch of the tree {Mh ~ W^'^Mq). Neverthe- 
less, the probability to have at least one visible progenitor 
in any merger history is very high (P = 68%). Interestingly, 
we found that the LAE candidates can be also observed 
as dropout galaxies since their UV magnitudes are always 
Muv < -18. 

On average the identified LAE stars have intermediate 
ages, t, « 150 — 400 Myr, and metallicities Z « 0.3-1^0; an 
exception is represented by five newly formed galaxies, which 
are hosted by small DM halos, Mh ~ W^Mq, which are yet 



visible as (faint) LAEs (La 
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their high star formation rate and extremely young stellar 
population, f» < 5 Myr. The low metallicity of these young 
galaxies, Z ^ 0.016 - 0.044Zq, reflects that of the MW 
environment at their formation epoch. Although rare (5 out 
of 80 LAEs), these Galactic building blocks could be even 
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more unambiguously identified among the least luminous 
LAEs, due to their larger Lya equivalent widths, EW = GO- 
ISO A with respect to those (f* 40 A) shown by older LAEsH . 
These small and recently virialized halos {zvir ~ 6) could be 
the progenitors of Fornax-like dwarf spheroidal galaxies (see 
Fig. 1 of Salvadori & Ferrara 2009). By identifying these 
LAEs, therefore, it would be possible to observe the most 
massive dSphs of the MW system just at the time of their 
birth. 

Uncertainties remain on the treatment of dust in calcu- 
lating the LF and observed properties of the MW progen- 
itor LAEs identified here, especially at the low luminosity 
end of the LF. Several aspects require additional study. As 
gas, metal and dust are preferentially lost from low mass 
halos, pushing the mass resolution of simulations to even 
lower masses would be important. Also, the amount of dust 
lost in SN-driven winds remains very poorly understood, an 
uncertainty that propagates in the evaluation of the con- 
tinuum and Lya radiation escaping from the galaxy. Fi- 
nally, Lya photons could be affected considerably by the 
level of dust clumping. It is unclear to what extent these ef- 
fects may impact the visibility of LAEs, as discussed in e.g. 
Dayal, Ferrara & Saro (2010). Progress on these issues is 
expected when high-resolution FIR/sub-mm observations of 
LAEs with ALMA will become available in the near future 
(Finkelstein et al. 2009; Dayal, Hirashita & Ferrara 2010). 
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